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Abstract. NML Cyg is a highly evolved OH/IR supergiant, one of the most prominent infrared objects due to its strong obscu- 
ration by dust, and supposed to be among the most luminous supergiants in the galaxy. We present the first diffraction-limited 
2.13 /im observations of NML Cyg with 73 mas resolution. The speckle interferograms were obtained with the 6 m telescope at 
the Special Astrophysical Observatory, and the image reconstruction is based on the bispectrum speckle-interferometry method. 
The visibility function declines towards the diffraction limit to ~ 0.6. 

Radiative transfer calculations have been carried out to model the spectral energy distribution, given by ground-based photo- 
metry and ISO spectroscopy, and our 2.13 pm visibility function. Additionally, mid-infrared visibility functions at 11 /im were 
considered. The observed dust shell properties do not appear to be in accordance with standard single-shell (uniform outflow) 
models but seem to require multiple components. Considering previous periods of enhanced mass-loss, various density en- 
hancements in the dust shell were taken into account. An extensive grid of models was calculated for different locations and 
strenghts of such superwind regions in the dust shell. To match the observations from the optical to the sub-mm domain requires 
at least two superwind regions embedded in the shell. The best model includes a dust shell with a temperature of 1000 K at 
its inner radius of 6.27?„, a close embedded superwind shell extending from 15.5R, to 21. 7R* with an amplitude (factor of 
density enhancement) of 10, and a far-out density enhancement at 186i?» with an amplitude of 5. The angular diameters of the 
central star and of the inner rim of the dust shell amount to 16.2 mas and 105 mas, resp. The diameter of the embedded close 
superwind region extends from 263 mas to 368 mas, and the inner boundary of the distant superwind region has a diameter of 
3'.' 15. In the near-infrared the dust condensation zone is limb-brightened leading to a corresponding ring-like intensity distri- 
bution. The grain sizes, a, were found to be in accordance with a standard distribution function, n(a) ~ a~ 3 ' 5 , with a ranging 
between o m i n = 0.005 /im and a max = 0. 15 ^m. The bolometric flux amounts to Fboi = 3.63 • Wm~ 2 corresponding to a 
central-star luminosity of L/Lq = 1.13 ■ 10 s ■ (d/kpc) 2 . Within the various parts of the dust shell, 1/r 2 density distributions 
could be maintained differing only in their amplitude A. A slight improvement of the far-infrared properties can be obtained if a 
shallower density distribution of p ~ 1/r 1,7 is considered in the distant superwind region. The present-day mass-loss rate was 
determined to be M = 1.2 • 10~ 4 M0/yr. The inner embedded superwind shell corresponds to a phase of enhanced mass-loss 
(with amplitude 10) in the immediate history of NML Cyg which began 59 yr ago and lasted for ~ 18 yr. Correspondingly, the 
outer superwind region is due to to a high mass-loss period (amplitude 5) which terminated 529 yr ago. 

Key words. Techniques: image processing - Circumstellar matter - Stars: individual: NMLCyg - Stars: mass-loss - Stars: 
supergiants - Infrared: stars 



1. Introduction 

The star NMLCyg (= V 1489 Cyg = IRC +40448) belongs 
to the most prominent infrared objects of the northern hemi- 
sphere. It was discovered by Neugebauer, Martz & Leighton 
(1965) being an extremely red object in their infrared survey. 
This redness is due to NMLCyg's strong obscuration by cir- 
cumstellar dust (e.g. Herbig & Zappala 1970[ Hyland et al. 
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1972). The spectra l type has been determined to be M6III 
(Herbig & Zappala 1970 ), however als o lum inosity cla sses o f 
up to la have been derived (e.g. Johnson 1968 Low et al. 1970 ). 



Various distance estimates exist for NMLCyg ranging 
from 200 pc (Herbig & Za ppala |1970| ) up to 3400 pc (Rowan- 
Robinson & Harris 1983 ). The most likely distance appears 
to be ~ 2kpc due to the imm ediate proximity to the Cyg 
OB2 association (Morris & Jura |1983[ ). N ML Cyg is partially 
surrounded by an HII region (Goss et al. 1974, Gregory & 
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Seaquist 1976, Habing et al. 1982) which is supposed to be 
the result of ionizing radiation from the Cyg OB2 association 
(Morris & Jura (19831). 



NMLCyg is associated with OH (Wilson et al. 
Masheder et al. 



1974, Benson & Mutel 1979, Bowers et al 
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1983), SiO (Snyder & Buhl |1975| , Dickins on et a l. |1978| , Morris 
et al. |1979[ ) and H 2 (Schw artz & Barrett |1970[ Schwartz et al. 



ity for the circumstellar shell w as me asured to be 21 km/s in 
the CO(2-l) line (Knapp et al. |l982|), 24 km/s in H 2 emis- 



1986), 26 km/s in thermal 



1992), and 28 km/s in OH 



sion profiles at 22 GHz (Engels et al. 
SiO emission at 87 GHz (Lucas et al. 
emission profiles at 1612MHz (Bowers et al. |1983j). Mass-loss 
rates, based on the CO and OH measurements and scaled to a 
distance of 1.8 kpc, are 1.5 • lO~ 4 M /yr and 1.6 • lO~ 4 M /yr, 
resp. 

Thus NMLCyg can be believed to be a highly evolved 
OH/IR supergiant of very large luminosity which suffers from 
stron g mas s-loss and is highly enshrouded by dust. Morris & 
Jura (1983 ) estimated its luminosity to be 5 • 10 5 L at the dis- 
tance of 2 kpc (correspondig to an initial mass of ~ 40M Q ). 
Consequently NML Cyg appears to be among the most lumi- 
nous supergiants in the Galaxy (Richards et al. 
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There are indications that the dust shell of NML Cyg may 
possibly deviate from spherical symmetry. Strong polarization 
was found by, e.g., Forbes (1967), Kruzewski (1971 ), and Dyck 
et al. (1971). The dust shell morphology appears to be complex. 
For instance, Diamond et al. (1984) found the OH envelope to 
be asymmetric, and Richards et al. (1996) concluded the exis- 
tence of a bipolar outflow from interferometric observations of 
the H 2 maser. 

Previous infrared speckle inter ferom etric observatio ns hav e 
been reported by Si bille et al. ( |l979|) Dyck et al. ( |l984| ), 
Ridgway et al. ( |l986| ), Fix & Cobb ( |l988| ) and Dyck & Benson 
(1992). Recent results from mid-infrared i nterfer ometry are 
given in Mo nnier et al. ( |1997| ), Danchi et al. ( |1999| ), and Sudol 
et al. ( |l999t )- 

In this paper we present diffraction-limited 73 mas bispec- 
trum speckle-interferometry observations of the dust shell of 
NML Cyg. Radiative transfer calculations have been performed 
to model the spectral energy distribution as well as visibilities 
at different wavelengths. 

2. Near-infrared observations and data reduction 

The NML Cyg speckle interferograms were obtained with the 
Russian 6 m telescope at the Special Astrophysical Observatory 
on June 13 and 14, 1998. The speckle data were recorded with 
our NICMOS-3 speckle camera (HgCdTe array, 256 2 pixels, 
sensitivity from 1 to 2.5 /im, frame rate 2 frames/s) through 
a continuum filter with a centre wavelength of 2.13/im and 
a bandwidth of 0.02 /im. Speckle interferograms of the unre- 
solved star HIP 102098 were taken for the compensation of the 
speckle interferometry transfer function. The observational pa- 
rameters were as follows: exposure time/frame 25 ms; number 
of frames 5000 (2000 of NML Cyg and 3000 of HIP 102098); 
2.13 fim seeing (FWHM) ~l'/0; field of view 7'.'8x7'.'8; pixel 
size 30.52 mas. A diffraction-limited image of NMLCyg with 
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Fig. 1. Left: Two-dimensional 2.13 /an visibility function of 
NMLCyg shown up to the diffraction limit (see right panel). 
The bright central structure shows that the central object is sur- 
rounded by a resolved dust shell. Right: Azimuthally averaged 
2.13 /im visibility of NML Cyg with error bars for selected fre- 
quencies. 



Fig. 2. Contour plot 
of the bispectrum 
speckle interferometry 
2.13/im-image of 
NMLCyg. Contour 
levels are plotted for 
2%, 3%, 4%, and 5% 
of peak intensity, and 
from 7% to 97% of 
peak intensity in steps 
of 3%. The dashed 
contours (2%, 3%) 
refer to the noise floor. 
North is up and east is 
to the left. 
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Fig. 3. Azimuthally 
averaged radial plots 
of the reconstructed 
diffraction-limited 
2.13/im-images of 
NMLCyg (solid line) 
and of the unresolved 
star HIP 102098 



200 (dashed line). 



73 mas resolution was reconstructed from the speckle interfer- 
ograms using the bispectrum speckle-interferometry method 
(Weigelt |1977|, Lohmann et al. |1983|, Hofmann & Weigelt 



1986 



Weigelt 



1991). The bispectrum of each frame consisted 



of 



37 million elements. The modulus of the object Fourier 
transform (visibility) was determined with the speckle interfer- 
ometry method (Labeyrie 197C). 

Figure [j] shows the reconstructed 2.13 /im visibility func- 
tion of NMLCyg, whereas Fig. || presents a contour plot of 
the reconstructed image. There is only little evidence for devi- 
ations from spherical symmetry. The object appears elongated 



T. Blocker et al.: Bispectrum speckle interferometry observations and radiative transfer models of NMLCyg 



3 



at the 4% intensity level, which, however, is already close to 
the underlying noise floor, with a position angle of ~ 140°, 
However, it is noteworthy that this position angle is close to 
the one found by 5 /mi interferometric measurements (~ 140°- 
150°; McCarthy 1979) and to those found for the symme- 
try axes of the SiO (~ 140°; Boboltz & Marvel 2000), H 2 
(~ 130°; Richards et al. 1996), and OH (~ 150°; Masheder et 
al. 1974) maser distributions. The GauB fit FWHM diameter of 
the dust shell was determined to be 121 mas. Fig. ^displays the 
azimuthally averaged diffraction-limited images of NML Cyg 
and the unresolved star HIP 102098. 



3. Mid-infrared visibilities 

Whereas the if -band visibility traces the hot dust in the imme- 
diate vicinity of the star, visibility functions in the mid-infrared 
can give spatial information on the colder dust located further 
out in the circumstellar shell where the 9.7 /an silicate feature 
forms. In the instance of NMLCyg and its heavy obscuration 
by dust, a considerable fraction of the observed flux is radi- 
ated at around 9.7 /im and due to the thickness of the shell the 
silicate feature appears in absorption (see Fig. |B|). Accordingly 
complementary TV-band (11.2/tm) interferometry can supply 
important constraints for the astrophysical interpretation of the 
object by, e.g., dust-shell modelling. Fig. |] refers to corre- 
sponding observations in this regime and shows the visibil- 
ity functions at 11.2/xm (Fix & Cobb 1988), at 10.4/xm and 
11.4/tm (Dyck & Benson 1992), at 11.15/xm (Monnier et al. 
1997, Danchi et al. 1999), and at 11.5/tm (Sudol et al. 1999). 
In the following sections we will focus on the most recent ones, 
i.e. the highest resolution observation of Monnier et al. (1997) 
and Danchi et al. (1999) with the ISI heterodyne interfero- 
meter in 1993, 1994, 1996, and 1998 and the latest low spa- 
tial frequency measurement of Sudol et al. (1999) at the 2.3 m 
Wyoming Infrared Observatory in 1997. 



4. Spectral energy distribution 

Photometric and spectrometric data in the optical and infrared 
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Fig. 4. Mid-infrared visibility functions of NMLCyg at 
11.2/im (Fix & Cobb 1988), at 10.4/xm and 11.4/tm (Dyck & 
Benson 1992), at 11.15/tm (Monnier et al. 1997; Danchi et al. 
1999), and at 11.5/tm (Sudol et al. 1999) 



for NMLCyg are given, e.g., by Johnson et al. (1965), Johnson 



qi967p, Low et al J1970J), Hackwell (|1972), Dycke t al. (|1974J), 
Strecker & Ney ( |1974| ) and Merrill & Stein ( |l976| ). NMLCyg 
shows small long-ter m var iations in its sp ectral energy distri- 
bution (Harvey et al. 1974 , Strecker 1975 ) with a period of ~ 
1000 days. Recently, Monnier et al. ( 1997 ) presented long-term 
10.2 /tm photometry covering a time basis of more than 12yr. 
They determined an amplitude of ~ 0.5 mag and a mean period 
of ~ 940 days. Thus NML C yg can be co nsidered as a semireg- 
ular variable. Monnier et al. ( 1997 , 1998 ) investigated the tem- 
poral variation of the mid-infrared spectra by means of 8- 1 3 /im 
spectrophotometry from 1991 to 1995 finding an almost con- 
stant spectral shape in this wavelength regime. Submillimeter 
(400 /im) an d milli meter (1.3 mm) obser vation s are given by 
Sopka et al. ( |1985| ) and Walmsley et al. ( |1991| ), resp. Besides 
these ground-based observations NML Cyg has also been ob- 
served by IRAS (Infrared Astronomical Satellite) in 1973 and 
by ISO (Infrared Space Observatory) in 1996. 

NML Cyg is highly reddened by t he inte rstellar medium 
and the circumstellar shell. Low et al. ( 1970 ) estimated a to- 
tal extinction of Ay tal « 9™5 for NML Cyg. Based on radio 
continuum observations Gregory & Seaquist (1976) obtained 
an interstellar extinction of Ay ea 4™6 whereas Lee (1970) es- 



timated E(B — V) = l'. n 2 for the interstellar contribution. We 
will use the latter value of E{B - V) = 1™2. This interstellar 
reddening was taken into account by adopting the method of 
Savage & Mathis ( |1979[ ) with A v = 3AE(B - V). 

The spectral energy distrubion (SED) of NMLCyg is 
shown in Fig. || and exhibits a 9.7 /im silicate profile in ab- 
sorption. It includes the (dereddened) data of Johnson et al. 
( 1 965), Dyck et al. ( 1 974) and Merrill & Stein ( 1 976) as well as 
the UKIRT spectrophotometry of Monnier et al. (1997). Note 
that the Merrill & Stein (1976) data has been scaled to match 
the UKIRT spectrophotometry as in Monnier et al. (1997), 
The spectral shapes of both data sets agree well. The lower 
panel of Fig. ^ shows the ISO-SWS spectroscopy (Justtanont 
et al. 1996) in comparison with the UKIRT data and the IRAS 
photometry. The ISO data was adjusted to match the ground- 
based photometry (see also Justtanont et al. 1996). In the long- 
wavelength region, the IRAS fluxes seem to be somewhat 
smaller than those provided by ISO. 



5. Dust shell models 

In order to model both the observed SED and the 2.13 /tm and 
11.15 /im visibilities, radiative transfer calculations for the dust 
shell were conducted. For the robust and non-ambiguous con- 
struction of dust shell models it is essential to take diverse and 
independent observational constraints into account. Apart from 
matching the spectral energy distribution, the consideration of 
spatially resolved information plays a crucial role for obtaining 
a reliable model (see e.g. Blocker 1999). 
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Fig. 6. Optical depth dependence of 
model SED (top) and model visibilities 
at 2.13 [im (bottom left) and 11.15 /im 
(bottom right) for T cff = 2500 K 
and T x = 1000 K. The optical depth 
r(0.55/im) ranges between 10 and 90. 
The calculations are based on a black 
body, Ossenkopf et al. (1992) silicates, 
a Mathis et al. (1977) grain size distri- 
bution with a max = 0.25, and a 1/r 2 
density distribution. Such a standard 
model does not appear to applicable to 
NML Cyg since different optical depths 
are required to match the various obser- 
vations. The symbols refer to the ob- 
servations (see text, Figs, [l], [| and ||) 
corrected for interstellar extinction of 
A v = 3.7: 

SED: *: Johnson et al. (1965), A: Dyck 
et al. (1974), o: Merril & Stein (1976), 
thick line: Justtanont et al. (1996) [ISO 
data]. The inlet shows the ISO data of 
Justtanont et al. (1996) (thick line) and 
the spectrophotometry of Monnier et al. 
(1997) (circles). 

Visibilities: +: 2.13 /im (present pa- 
per); 

A: 10.4 ^tm (Dyck & Benson 1992), 
T: 11.4/xm (Dyck & Benson 1994), 
o: 11.2 /an (Fix& Cobb 1988), *: 
11.15 /im (Monnier et al. (1997), v : 
11.15/im Danchi et al. (1999), □: 
11.5 Aim Sudol et al. (1999). 



5. 1. The radiative transfer code 



We used the code DUSTY developed by Ivezic et al. (1997), „ 



which solves the radiative transfer problem in spherical sym- 
metry. It considers absorption, emission and scattering and uti- 
lizes the self-similarity and scaling behaviour of IR emission 
from radiatively heated dust (Ivezic & Elitzur 1997 ). To solve 
the radiative transfer problem the spectral shape of the central 
source's radiation has to be specified and various properties of 
the surrounding dust shell are required, viz. (i) the chemical 
composition and grain size distribution; (ii) the dust tempera- 
ture at the inner boundary; (iii) the relative thickness, i.e. the 
ratio of outer to inner shell radius; (iv) the density distribution; 
and (v) the total optical depth at a given reference wavelength. 
The code has been expanded for the calculation of synthetic 
visibilities as described by Gauger et al. (1999). 



5.2. Single-shell models 

Dust shell models for NMLCyg were reported in Rowan- 
Robinson & Harris (1983), Rigdway et al. (1986), and Monnier 



et al. (1997). The first two studies used uniform outflow mod- 
els, i.e. approximated the density distribution of the dust shell 
r~ 2 . Monnier et al. (1997) found that uniform outflow 
models appear not to be consistent with their near- and mid- 
infrared observations and proposed a double-shell structure for 
the dust envelope. It should be kept in mind, however, that the 
older calculations are partly based on different input quantities, 
as optical constants, grain sizes, or shape of the central source 
radiation. 

We calculated various uniform outflow models considering 
the following parameters within the radiative transfer calcula- 
tions: SED and visibilities were modelled for effective tem- 
peratures of T c ff = 2000 to 3000 K and black bodies as cen- 
tral sources of radiation. We used the silicates of Ossenkopf 
et al. (1992) which were found to be well suited to fit the sil- 



icate feature by Monnier et al. (1997). Some test calculations 
were p erformed with the optical const ants of of Draine & Lee 
( 1984 ) and David & Pegourie ( 1995 ) as well. For the grain- 
size distribution we consid ered grain size distributions accord- 
ing to Mathis et al. ( |1977| , hereafter MRN), i.e. n(a) ~ a~ 3 5 , 
with 0.005 fim < a < (0.13 to 0.45) /im as well as single- 
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Fig. 5. Spectral energy distribution of NML Cyg (corrected for 
interstellar extinction of A v =3.7). Top: Photometry of Johnson 
et al. (1965), Dyck et al. (1974) and Merrill & Stein (1976) as 
well as the UKIRT spectrophotometry of Monnier et al. (1997). 
The Merrill & Stein (1976) data was scaled to match the 
UKIRT photometry. Bottom: ISO-S WS spectrum of Justtanont 
et al. (1996), UKIRT spectrophotometry of Monnier et al. 
(1997), and IRAS photometry. The ISO data was adjusted to 
match the ground-based photometry (see also Justtanont et al. 
1996). 

sized grains with a — 0.05 to 0.3 /im. The shell thickness 
^jut = r out /ri amounted to 10 3 to 10 4 with r out and t\ be- 
ing the outer and inner radius of the shell, respectively. With 
these quantities specified, the remaining fit parameters are the 
dust temperature at the inner boundary, Ti, which, in turn, de- 
termines the radius of the shell's inner boundary, t\, and the 
optical depth, r, at a given reference wavelength, A re f . We refer 
to A re f = 0.55 /im. Models were calculated for dust tempera- 
tures between 400 and 1500 K and optical depths between 10 
and 100. 

The neccessity to consider various observational con- 
straints in order to obtain dust-shell properties by radiative 
transfer calculations is demonstrated in Fig. |[ It shows that 
the SED and the near-infrared and mid-infrared visibilities 
cannot be fitted simultaneously by a given set of parameters. 
Though each single observation can be matched reasonably 



well, different values for the optical depth are required. For ex- 
ample, the K-band visibility requires r(0.55/im) = 10, the 
N-band visibility T(0.55/im) = 30, and the silicate feature 
r(0.55/im) = 50 for the chosen dust-shell properties. High 
optical depths (r(0.55/xm) ^ 30) lead to a flux deficit in the 
optical. 

This example refers to a given dust temperature of T\ = 
1000 K at the inner boundary and an MRN grain-size distri- 
bution n(a) with a max = 0.25 /im. Additionally the mod- 
els rely on the assumption of a standard density distribution 
p{r) ~ 1 /r 2 whose validity has been questioned by Monnier 
et al. (1997). Before abondaning the uniform outflow model, 
however, one has to consider the effects of changes of the for- 
mer two dust-shell properties, T\ and n(a). 

The if-visibility is very sensitive against scattering and, 
thus, depends strongly on the assumed grain sizes. Increasing 
the (maximum) grain size has a similar effect as increasing 
the optical depth and results into a stronger decline of the K- 
visibility (cf. e.g. Blocker et al. 1999). This holds both for 
single-sized grains and grain-size distributions. If one increases 
the maximum grain-size of an MRN-like distribution most par- 
ticles are still small due to the steep decline of the distribution 
function towards larger grain sizes. However, due to the very 
strong dependence of the absorption and scattering properties 
on the grain sizes, these "few" larger particles determine the 
if -visibility equally as much as the bulk of the smaller grains 
(Kriiger & Sedlmayr 1997, Winters et al. 1997). On the other 
hand, at longer wavelengths the changes induced by a larger 
maximum grain size are predominantly of minor nature. For in- 
stance, the SED silicate feature and the mid-infrared visibility 
are only mildly affected. This fact allows, in principle, to fix the 
optical depth by matching the silicate feature and mid-infrared 
visibilities (if possible) and to fix independently the grain-size 
distribution by the near-infrared data. The impact of maximum 
grain sizes of MRN-like distribution on the various observa- 
tional quantities is illustrated in Fig. 0. These calculations show 
that the low-frequency part of the if-band visibility may be 
fitted for smaller optical depths if a max is sufficiently small. 
However, the silicate feature and the SED's long-wavelength 
tail are then only poorly matched. The most favorable maxi- 
mum grain size within an MRN-like distribution was found to 
be close to a max = 0.15 fj,m. 

To complete the discussion on grain sizes it has to be con- 
sidered how an MRN-like distribution compares with single- 
sized grains and which choice is better suited for the instance 
of NMLCyg. Whereas the A"-visibility behaves similarly in 
both cases, i.e. strongly declines with for increasing grain sizes, 
SED and mid-infrared visibilities are more affected by varia- 
tions of single-sized grains than by those of a distribution. The 
interplay between optical depth and grain size is now much 
stronger. In summary, we found that single-sized grains appear 
to give worse overall fits compared to grain-size distributions. 

Finally, we want to address the question of the dust tem- 
perature at the inner shell boundary. The hotter the dust at the 
inner boundary of the shell, the more flux is radiated in the 
near-infrared and the less dominant is the long-wavelength tail 
of the SED. Concomitantly, the visibilities decline at low spa- 
tial frequencies the steeper, the hotter the inner dust-shell rim. 
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Fig. 7. Maximum grain-size depen- 
dence of model SED (top) and model 
visibilities at 2.13 /im (bottom left) and 
11.15 /im (bottom right) for T c g = 
2500 K, Ti = 1000 K, r(0.55^m) = 
30. The maximum grain size a max 
ranges between 0.15 /im and 0.45 [im. 
The calculations are based on a black 
body, Ossenkopf et al. (1992) silicates, 
a Mathis et al. (1977) grain size distri- 
bution, n(a) oc a~ 3 - 5 , and a \/r 2 den- 
sity distribution. The symbols refer to 
the observations (cf. Figs, [j], Q and ^) 
corrected for interstellar extinction of 
A v = 3.7 and are the same as those in 
Fig- 1 



Inspection of the model grid shows that the dust temperature 
at the inner boundary should be close to 1000K. Somewhat 
higher values, e.g. 1200 K, are possible as well. 

All parameter variations discussed so far (of optical depths, 
dust temperatures, grain sizes etc.) failed to give a consistent 
overall fit, thus the only basic parameter left to vary is the den- 
sity distribution (assuming spherical symmetry). Consequently, 
it can be concluded that the assumption of a uniform out- 
flow model does not hold for NML Cyg and that, for instance, 
steeper density distributions or multiple shells have to be con- 
sidered. A similar conclusion was drawn by Monnier et al. 
(1997) who introduced two shells with exponentially decreas- 
ing density distributions in order to match the observations. 

5.3. Multiple dust-shell components 

The evolution of massive supergiants can be accompanied 
by stages of episodic mass-loss events as superwind phases 
or even mass-loss outbursts, as for example, during the Red 
Supergiant Branch or the Luminous Blue Variable stage (see, 
e..g., Vanbeveren et al. 1994, Heger et al. 1997, de Jager & 
Nieuwenhuijzen 1997, Langer et al. 1998). During transitory 
phases of heavy mass-loss the density in the dust shell is in- 
creased. This leads to regions in the dust shell with density en- 
hancements over the normal 1 /r 2 distribution. The radial den- 
sity distribution may also change within such superwind shells. 



For more details, see Suh & Jones (1997). For example, such 
superwind models have successfully been applied to the hy- 
pergiant IRC + 10 420 that is believed to have already moved 
off the RSG currently evolving into the Wolf-Rayet stage (see 
Blocker et al. 1999). 

Concerning the introduction of superwind models different 
degrees of complexity can be considered: (i) single jumps with 
enhancement factors, or amplitudes, A at radii Y n i = r n \/r\; 
(ii) a single superwind shell with amplitude A ranging from Y n i 
to Y n 2 \ (iii) various shells; and (iv) jumps or shells as in (i)-(iii) 
but with density distributions different from the 1/r 2 uniform 
outflow model. Examples of the superwind models (i)-(iii) are 
illustrated in Fig. ^. 

In the following sections such superwind models will be 
discussed in more detail. In each case the suitability of other 
dust-shell parameters was carefully checked. We found a op- 
tical depth of r(0.55/im)=30 to fit the silicate feature. Within 
an MRN-like grain-size distribution a maximum grain-size of 
Omax = 0.15/im was chosen in order to match the overall pat- 
tern of the SED and the low-frequency curvature of the X-band 
visibility. We stayed with the optical constants of Ossenkopf et 
al. (1992) due to their excellent match of the silicate feature's 
shape. The dust-temperature at the inner dust-shell boundary 
amounted to 1000K, the effective temperature of the central 
star to 2500 K. 
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Fig. 8. Relative density distribution for different superwind 
models with p ~ 1/r 2 . The upper panel (a) shows a density 
enhancement at Y n i = r n x/fi — 2.5, i.e. nearby the inner 
dust-shell rim, with an amplitude of A = 10. The middle panel 
(b) refers to a superwind shell ranging from Y n \ — 2.5 to 
Y n 2 = 3.5 with A = 10, and the bottom panel (c) shows this 
superwind shell and a additional /ar-out density enhancement 
at Y f i = 30 with A = 5. 

5.3.1 . Two-component shells 

The most simple model of a previous epoch of enhanced mass- 
loss is that of a single density enhancement as shown in the 
upper panel of Fig. ||. Calculations were conducted for rel- 
ative distances Y n i between 1.5 and 6.5 and amplitudes A 
ranging from 10 to 20. The best model found is that with 
Y n i = 2.5 and A — 10 which fits the SED from the opti- 
cal to the infrared and matches the if -band visibility reason- 
ably well. The inner boundary of the dust shell amounts to 
r\ = 5.8i?* (Ti = 1000 K) corresponding to an angular diam- 
eter of 95 mas. At Y n \ = 2.5 the dust temperature has dropped 
to 621 K. More distant density jumps lead to worse fits of the 
silicate feature and the far-infrared regime, less distant ones 
are not in agreement with the Jf-band visibility (for given am- 
plitude). However, the mid-infrared visibility was only poorly 
matched by all of these models giving too extended shells at 
this wavelength. This, in turn, indicates that the cold dust-shell 
regions which determine the mid-infrared visibility have too 
high densities within this model. 

5.3.2. Embedded superwind shells 

Since the if -band visibility traces the hot dust, it is sensi- 
tive to the place and strength of the density jump but reacts 
only moderately to the extension of the density enhancement. 
Consequently, the above model can be improved by introducing 
an embedded high-density or superwind shell of given exten- 



sion (see Fig. || middle panel) instead of a single jump. Since 
the inner boundary and the amplitude of such a shell is strongly 
constrained by the K -band visibility, the inner boundary was 
fixed to Y n i = 2.5 and the amplitude to A = 10 within the cor- 
responding model grid. Additional test calculations with other 
values for Y n i and A were performed as well. The outer bound- 
ary Y n 2 was varied between 3.0 and 10.5. SED and visibility 
models for different values of Y n 2 are shown in Fig. ^. It turned 
out that only a close and comparatively thin shell ranging from 
Yni = 2.5 to Yni = 3.5 appears to be in agreement with 
both the visibility data at different wavelengths and the SED 
from the optical to the mid-infrared. Now, the inner boundary 
of the dust shell is located somewhat more outwards, viz. at 
Ri = 6.2i?* corresponding to an angular diameter of 101 mas. 
The superwind shell diameter extends from 252 to 353 mas, the 
corresponding dust temperature at the inner and outer bound- 
ary amount to 587 K to 41 1 K, resp. However, as obvious from 
Fig. ^ the model appears to be still incomplete due to a lack of 
sufficiently high fluxes at longer wavelengths (A <; 20^m). 

5.3.3. Two superwind shells 

The above model of a thin density-enhanced region close to 
the inner boundary of the dust shell gives a reasonable fit to 
the observed properties from the optical to the mid-infrared. 
However, due to the small extension of the embedded super- 
wind shell, the densities in the outer regions of the dust shell 
appear to be too small to provide sufficiently high far-infrared 
fluxes. This may indicate either the existence of a second, far- 
out density enhancement or a shallower radial density depen- 
dence (e.g. p ~ 1/r 15 ). In the case of a second episodic mass- 
loss phase, the corresponding density enhancement should be 
at moderate distances, close enough to raise the fluxes at A > 
20pm but distant enough in order to preserve the properties of 
the radiative transfer models at shorter wavelengths. 

Grids were calculated for such an additional far-out super- 
wind shells with an inner boundary located between Yfi = 7 to 
50 and amplitudes A ranging between 2 and 50. The best model 
found was that with Yfi = 30 and a moderate amplitude of 
A = 5. For simplicity we considered here only density jumps 
as shown in Fig. [| Note that though the existence of such 
a distant density enhancement appears to be well constrained 
by the observations within the current model, the correspond- 
ing properties of the distant shell are less constrained. For in- 
stance, different combinations of location, amplitude and den- 
sity slope can lead to models of comparable quality. Hovever, 
shells at significantly smaller or larger distances can clearly be 
excluded. 

Since the introduction of a far-out density enhancement 
hardly affects the inner regions, the inner boundary of the dust 
shell keeps at 6.2 stellar radii corresponding to an angular 
diameter of 105 mas. The angular diameter increased some- 
what compared to the case without an outer shell due to the 
the (slight) increase of the bolometric flux. The diameter of 
the embedded close superwind shell extends from 263 mas to 
368 mas, and the inner boundary of the distant shell has a diam- 
eter of 3'.' 15. The dust temperature within the embedded super- 
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Fig. 9. Model SED (top, left) and 
model visibilities at 2.13 /itn (bot- 
tom left) and 11.15 pm (bottom right) 
for T cf f = 2500 K, r(0.55/im)=30, 
Ti=1000K, and different embedded 
superwind shells with amplitude ^4=10 
extending from Y n i = 2.5 to Y n 2 = 
3.0-8.5. The density distribution is ex- 
emplified for "5^1=2.5 and 1^2=3.5 in 
the above panel (top, right) (cf. Fig. 
The calculations are based on a black 
body, Ossenkopf et al. (1992) silicates, 
and a Mathis et al. (1977) grain size 
distribution with a max =0.15/im. The 
symbols refer to the observations (cf. 
Fig.[l|,|and|) corrected for interstellar 
extinction of A v =3.7 and are the same 
as those in Fig. 0. 



wind shell ranges from 584 K to 430 K and within the distant 
shell from 149 K to 40 K. 

Fig. |ll| gives the fractional flux contributions of the emerg- 
ing stellar radiation, of the scattered radiation and of the ther- 
mal dust emission as a function of wavelength for this model. 
At 2.2/im, the flux is mainly determined by direct stellar light 
(67%). Scattered radiation and dust emission contribute 6% 
and 27%, resp. Thus, roughly 80% of the circumstellar shell's 
iT-band radiation is due to thermal dust emission. At wave- 
lengths larger than 10pm dust emission completely dominates. 
The bolometric flux amounts to F^ \ = 3.63 • 10~ 9 Wm~ 2 . 
Accordingly, the central star has a luminosity of L/Lq = 
1.13 • 10 5 • (d/kpc) 2 and an angular diameter of 9* = 1.74 • 
10 9 y/F hol /T? s ~ 16.2 mas. 

Fig. 12 shows the normalized intensity distribution at 



2.13/im (upper panel) and 11.15/im (lower panel) as func- 
tion of angular distance. The barely resolved central peak cor- 
responds to the central star (angular diameter: 16.2 mas). At 
2.13/im the inner boundary of the dust shell (radius: 52.5 mas) 
is limb-brightened resulting in a ring-like intensity distribution. 
The optical depth (r(2.13^m) = 1.5) is still moderate lead- 
ing to a noticable effect of a limb-brightened dust condensation 
zone (cf. Ivezic & Elitzur 1996, Blocker et al. 1999). In the N 
band (r(11.5/xm) = 2.2) this is hardly visible. The embedded 
superwind shell leads to an increase of the intensity distribu- 
tion's waist at angular radii between 131.5 mas and 184 mas. 



Finally, the far-out shell provides broad low-intensity wings at 
angular radii larger than ~ l'/5. 

Asuming an outflow velocity of 25km/s and a distance 
of 1.8 kpc, the kinematical ages of these shells can be cal- 
culated. The inner embedded shell corresponds to a phase of 
enhanced mass-loss lasting for ~ 18 yr in the immediate his- 
tory of NMLCyg. It began 59.2 yr ago and ceased 41.5 yr ago. 
During this time span the mass-loss rate increased by a factor 
of 10. Correspondingly, the outer shell, i.e. the density jump 
at Yfi = 30, is due to to a high mass-loss period which ter- 
minated 529.3 yr ago. Throughout this phase in the past matter 
was expelled with a fivefold increased rate. With a dust-to-gas 
ratio of 0.005 and a specific dust density of 3 g cm~ 3 the model 
gives a present-day mass-loss rate of M = 1.2 • 10~ 4 M Q /yr 
which is in line with the CO and OH observations (see Sect. [[]). 

5.3.4. Other density slopes 

In the above sections the uniform-outflow model was aban- 
doned due to its mismatch of the observed properties. The in- 
troduction of two superwind shells, or density enhancements 
of certain amplitude A, leads to a model which fits the ob- 
servations from the optical to the infrared. However, so far it 
was assumed that p ~ 1/r 2 holds within each dust-shell re- 
gion. Assuming the outflow velocity kept constant any shal- 
lower (steeper) density slope than p ~ 1/r 2 means that the 
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Fig. 10. Model SED (top) and model 
visibilities at 2.13 /mi (bottom left) and 
11.15 /mi (bottom right) for an inner 
superwind shell of amplitude ^4=10 be- 
tween y n i=2.5 and 3^2=3.5, and an ad- 
ditional far-out density enhancement at 
Yfi=30 with A=5. The density distri- 
bution is displayed in the above panel 
(top, right) (cf. Fig. fy. The models 
refer to T cff =2500K, Ti=1000K and 
optical depths t(0.55/xto) ranging be- 
tween 10 and 90. The calculations are 
based on a black body, Ossenkopf et 
al. (1992) silicates, and a Mathis et 
al. (1977) grain size distribution with 
a max =0.15/im. The symbols refer to 
the observations (cf. Fig. [FJ, |] and ^) 
corrected for interstellar extinction of 
A v =3.7 and are the same as those in 
Fig- 1 



mass-loss rate has decreased (increased) with time in the re- 
spective dust-shell regions. 

The double-shell superwind models still seem to predict a 
somewhat larger spectral index, i.e. faster flux decrease, than 
the ISO observations indicate for A <; 20 (Fig. |l0|), and than 
400/mi and 1.3 mm data of Sopka et al. (1985) and Walmsley et 
al. (1991), resp., show. Additional models with shallower den- 
sity distributions in the outer shell were calculated to improve 
the model. A somewhat shallower distribution of p ~ l/r 1 - 7 
in the far-out regions (Y > 30) of the dust shell appears to be 
slightly better suited than the standard model with p ~ l/r 2 
and matches also the 400/im observations. The 1.3mm data, 
however, cannot be matched by any of the models. It has al- 
ready been noted by Sopka et al. (1985) and Walmsley et al. 
(1991) that the fluxes at 400/im and 1.3mm are in excess of 
standard models. It has been argued that flat density distribu- 
tions (as introduced above) or an increased emissivity of the 
silicates may be responsible for this excess. 

If one considers different density slopes in the outer shell 
to improve the model, the question arises if this should not be 
considered for the whole dust shell model in general. For in- 
stance, Monnier et al. (1999) found two shells with exponen- 
tially decreasing density distributions appropriate to match the 
observations. In order to prove (or disprove) the applicability of 
the assumed 1 /r 2 distribution in the inner regions of the shell 
(Y < 30), we re-calculated the main body of the grid presented 



so far including the superwind shells for density distribution 
ranging from p ~ l/r to p ~ l/r 7 and combinations of it. 
We found the same best model as presented above to match 
the observations from the optical to the sub-mm domain, i.e. 
a (close-to) 1 /r 2 overall density distribution with a thin super- 
wind shell (Y — 2.5 to 3.5) near to the inner dust-shell rim and 
a far-out density enhanced region (Y J; 30). Shallower distri- 
butions in the inner region can be excluded, for instance due 
to the mismatch of the mid-infrared visibility. Steeper distribu- 
tions partly match the observations but were found to be unable 
to provide an overall fit. If one considers only the embedded su- 
perwind shell, models can be constructed which show, in prin- 
ciple, similar properties as the one presented above but have a 
much steeper density distribution within the embedded super- 
wind shell, i.e. as steep as l/r 5-7 . Then the superwind shell 
has to be somewhat more extended and to exhibit a larger am- 
plitude (e.g. A — 20). However, preference can be given to the 
model showing superwind shells of enhanced density but pre- 
serving a l/r 2 distribution in the various parts of the dust shell 
due to a better overall match of the observations. 

6. Summary and discussion 

The first diffraction-limited 2.13 /im observations of NML Cyg 
with 73 mas resolution were presented. The speckle interfer- 
ograms were obtained with the 6 m telescope at the Special 
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Fig. 11. Fractional contributions of the emerging stellar radia- 
tion as well as of the scattered radiation and of the dust emis- 
sion to the total flux as a function of the wavelength for the 
model with for an inner superwind shell of amplitude A=10 
between Yni=2.5 and 1^2=3.5 and an additional far-out den- 
sity enhancement at lfi=30 with A=5 (see Fig. |l0|). Model 
parameters are: black body, T cS = 2500 K, T a = 1000 K, 
T"o.55/im = 30, Ossenkopf et al. (1992) silicates, and a Mathis 
et al. (1977) grain size distribution with a max = 0.15 /zm. 



Astrophysical Observatory, and the image reconstruction was 
based on the bispectrum speckle-interferometry method. In 
order to interpret these observations radiative transfer calcu- 
lations have been carried out. We modelled the spectral en- 
ergy distribution, in particular the ISO LWS spectroscopy 
(Justtanont et al. 1996), the 2.13 /mi visibility function, and 
mid-infrared visibility functions at 11 /im (e.g. Danchi et al. 
1999). Preference was given to those observations conducted 
in phase or close to the 2.13 /im observations in June 1998. 

The observed dust shell properties do not appear to be in ac- 
cordance with standard single-shell models (uniform outflow) 
but seem to require multiple components. Since single-shell 
models fail to match the observations, density enhancements 
in the dust shell were considered. Such density enhancements 
correspond to previous periods of enhanced mass-loss. The 
episodic and/or stochastic occurence of such superwind phases 
seems to be a typical feature of the late stages of massive star 
evolution. An extensive grid of models was calculated for dif- 
ferent locations and amplitudes of such superwind regions in 
the dust shell. It turned out that at least two superwind phases 
have to be taken into account in order to to match the observa- 
tions from the optical to the sub-mm domain. The best model 
is that of a dust shell with a temperature of 1000K at its inner 
radius of 6.2i?„. It includes a close embedded superwind shell 
extending from 15.5-R* to 21.7i?* with an amplitude of 10 and 
a far-out density enhancement at 186i?* with an amplitude of 5. 
The inner rim of the dust shell, i.e. the dust condensation zone, 
has a dust temperature of T\ = 1000 K and an angular diameter 
of 105 mas. At 2.13 /im this zone is limb-brightened leading to 
a ring-like intensity distribution. The diameter of the embedded 
close superwind region extends from 263 mas to 368 mas and 
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Fig. 12. Normalized intensity at 2.13 /im (top) and 11.15 /im 
(bottom) vs. angular displacement $ for a superwind model 
with an inner embedded superwind shell of amplitude A = 10, 
located between Y a i = 2.5 and Y n 2 — 3.5, and an additional 
far-out density enhancement at Yfi — 30 with A = 5. The cen- 
tral peak belongs to the central star. The inner hot rim of the cir- 
cumstellar shell has a radius of 52.5 mas, the inner embedded 
superwind shell extends from 13 1.5 mas to 184 mas, and the 
cool far-out component is located at a radius of 1500 mas. All 
loci correspond to local intensity maxima or plateaus. In the K- 
band the inner rim of the circumstellar shell is noticeably limb- 
brightened. The models refer to T cS = 2500 K, T x = 1000 K, 
and r(0.55/im) = 30. The calculations are based on a black 
body, Ossenkopf et al. (1992) silicates, and a Mathis et al. 
(1977) grain size distribution with a max = 0.15/im. 



is composed of relatively cool dust with temperatures ranging 
from 583 to 430 K. At the inner boundary of the distant super- 
wind region the dust temperature has dropped to 149 K corre- 
sponding to an inner diameter of 3'.' 15. 

Fig. [l3] shows the 2.13 /im visibility for baselines up to 
36 m and different assumptions on the superwind regions in 
comparison with a standard uniform outflow model. The cen- 
tral star, having a diameter of ~ 16 mas, becomes partially re- 
solved for baselines larger than 10m and is completely resolved 
at a baseline of 3 1 m. 

During the entire grid calculations the validity of other pa- 
rameters as, optical depths, dust temperatures, grain sizes, ef- 
fective temperatures, etc., was comprehensively checked. The 
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Fig. 13. Model visibilities at 2.13^m in the long-baseline inter- 
ferometry range up to baselines of 36m (6 m: g=13.4 arc sec -1 ) 
for different assumptions on the density distribution. The solid 
line gives the best fit to the observations and refers to a dust 
shell with an inner embedded superwind region with amplitude 
A = 10, located between Y n \ = 2.5 and Y n 2 = 3.5, and an ad- 
ditional far-out density enhancement at Yfi = 30 with A = 5. 
The other lines correspond to models without an outer den- 
sity enhancement and a superwind shell of equal properties, of 
larger amplitude, and larger extension, resp. A standard model 
without any density enhancement is given as well. The models 
refer to T cff = 2500 K, Tj = 1000 K, and r(0.55^m) = 30 
and are based on a black body, Ossenkopf et al. (1992) sili- 
cates, and a Mathis et al. (1977) grain size distribution with 
a max = 0.15/im. 

grain sizes, a, were found to be in accordance with a stan- 
dard distribution function, n(a)<~a -3 ' 5 , with a ranging be- 
tween a m in = 0.005 yum and a max = 0. 15 fim. The effective tem- 
perature was chosen to be 2500 K. The temperature of a gi- 
ant with spectral type M6III is - 3200 K (Perrin et al. 1998) 
whereas supergiants (Til) have systematically lower tempera- 
tures (Dyck et al. 1996, 1998). For an M6 supergiant the ef- 
fective temperature can be estimated to be close to 3000 K. 
However, (moderate) variations of the effective temperature af- 
fect the radiative transfer models only slightly. Thus, though a 
somewhat higher effective temperature appears to better com- 
ply with the spectral type, differences to the presented T e g = 
2500 K grid are only of minor nature and results and conclu- 
sions of the presented study remain unchanged. Corresponding 
test calculations have been performed for T cS = 3000 K and 
T e ff = 3500 K. The temperature at the inner rim of the dust 
shell was found to be more or less fixed at 1000K in order 
to be in accordance with the observations. However, slightly 
higher temperatures, e.g. 1 100 K, appear to be possible as well. 
The bolometric flux amounts to Fboi = 3.63 • 10 -9 Wm" 2 cor- 
responding to a central-star luminosity of L/L Q = 1.13 • 10 5 • 
(d/kpc) 2 . Accordingly, NMLCyg is a true massive supergiant 
with L ~ 4 • 10 5 L Q at the distance of the Cyg OB2 association. 

Though the density distribution was scaled by amplitudes 
A in various parts of of the dust shell, its general 1/r 2 depen- 
dence proved to be appropriate to model the observations and 
thus could be maintained. Alternative models were found if one 



considers a much steeper density distribution in the inner em- 
bedded superwind region. Then, however, the superwind shell 
has the same inner diameter but a larger amplitude and exten- 
sion. A slight improvement of the far-infrared properties can 
be obtained if a shallower density distribution of p ~ 1/r 1,7 
is considered in the distant superwind region. Then, assuming 
that the outflow velocity kept constant, the mass-loss rate in 
the outer component has decreased with time. The present-day 
mass-loss rate can be determined to be M = 1.2- 10~ 4 M Q /yr. 
The inner embedded superwind shell corresponds to a phase 
of enhanced mass-loss (M = 1.2 • lO~ 3 M /yr) in the imme- 
diate history of NMLCyg which began 59.2 yr ago and lasted 
for ~ 18 yr if d = 1.8 kpc and v = 25km/s. Correspondingly, 
the outer superwind region is due to to a high mass-loss period 
(M = 6.0 • 10~ 4 M /yr) which terminated 529.3 yr ago. 

The inner diameter of the dust shell of 105 mas is close 
to the one found by the standard uniform-outflow models of 
Rowan-Robinson & Harris (1983) and Ridgway et al. (1986) 
(80 mas and 90 mas, resp.) but only half as large as the one 
found by Monnier et al. (1997) (240mas). Note, however, 
that the repsective models differ in various aspects. On the 
one hand, the simultaneous match of SED and visibilities ap- 
pear to require density distributions different from the stan- 
dard uniform-outflow assumption as indicated by the models 
of Monnier et al. (1997) and the present models. On the other 
hand, Monnier et al. (1997) used a two-component model with 
exponentially declining density distributions that is different 
from the current superwind model which, in turn, shows single 
density jumps but preserves 1/r 2 density courses in the various 
dust shell parts. The present study matches the observations 
for a broad wavelength range, i.e. from ~ 0.7^m to 400^m. 
Another difference concerns the choice of other dust-shell pa- 
rameters, in particular the optical constants, the dust-shell tem- 
perature and the dust-grain sizes. Furthermore, it has to be kept 
in mind that the observational data often belong to different 
epochs and in the case of long-baseline interferometry on the 
used position angle as well. For instance, Sudol et al. (1999) 
found their low spatial-frequency 11/im data not to be consis- 
tent with the higher frequency data of Monnier et al. (1997) 
possibly pointing to an expanding double-shell structure (due 
to the different epochs) or to the existence of bipolar structures 
(due to the different position angles). They stress the need of 
coeval data covering several position angles to solve this issue. 
Danchi et al. (1999) presented follow-up observations of the 
intial measurements of Monnier et al. (1997) giving evidence 
of the dust-shell expansion at 11pm. We note that it was essen- 
tial for the current superwind models to consider the Danchi 
et al. (1999) data being coeval with our if -band measurements 
instead of taking only the former Monnier et al. observations. 
Finally, all models of NML Cyg available so far rely on the as- 
sumption of spherical symmetry whose applicablity still has 
to be proven. Richards et al. (1996), for instance, discussed 
the existence of a bipolar outflow from interferometric maser 
observations. Nature is certainly much more complex than the 
simple models given in the present study and great care has to 
be taken during model construction concerning its reliability 
and unambiguousness. Nevertheless, it can be concluded that 
the dust shell of NML Cyg consists of a multiple-component 
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structure with two regions attesting the ocurrence of previ- 
ous superwind phases, one in its immediate past, i.e. decades 
ago, and one hundreds of years ago. Thus, the red supergiant 
NML Cyg proves to be an important example of episodic mass- 
loss phases in the evolution of massive stars along the Red 
Supergiant Branch. 
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